Robust multigrid preconditioners for cell-centered finite volume discretization of the high-contrast diffusion equation
نویسندگان
چکیده
We study a conservative 5-point cell-centered finite volume discretization of the high-contrast diffusion equation. We aim to construct preconditioners that are robust with respect to the magnitude of the coefficient contrast and the mesh size simultaneously. For that, we prove and numerically demonstrate the robustness of the preconditioner proposed by Aksoylu et al. (2008, Comput. Vis. Sci. 11, pp. 319–331) by extending the devised singular perturbation analysis from linear finite element discretization to the above discretization. The singular perturbation analysis is more involved than that of finite element because all the subblocks in the discretization matrix depend on the diffusion coefficient. However, that dependence is eliminated asymptotically. This allows the same preconditioner to be utilized due to similar limiting behaviours of the submatrices; leading to a narrowing family of preconditioners that can be used for different discretizations—a desirable preconditioner design goal. We compare our numerical results to standard cell-centered multigrid and observe that performance of our preconditioner is independent of the utilized prolongation operators and smoothers. As a side result, we also prove that the solution over the highly-diffusive island becomes constant asymptotically. Integration of this qualitative understanding of the underlying PDE to our preconditioner is the main reason behind its superior performance. Diagonal scaling is probably the most basic preconditioner for highcontrast coefficients. Extending the matrix entry based B. Aksoylu (corresponding author) · Z. Yeter Department of Mathematics & Center for Computation and Technology, Louisiana State University 216 Johnston Hall, Baton Rouge LA, 70803 USA E-mail: [email protected], [email protected] spectral analysis introduced by Graham and Hagger, we rigorously show that the number of small eigenvalues of the diagonally scaled matrix depends on the number of isolated islands comprising the highly-diffusive region.
منابع مشابه
Robust preconditioners for the high-contrast Stokes problem
We study the Stokes equation with high-contrast viscosity coefficient and this regime corresponds to a small Reynolds number regime because viscosity is inversely proportional to the Reynolds number. Numerical solution to the Stokes flow problems especially with high-contrast variations in viscosity is critically needed in the computational geodynamics community. One of the main applications of...
متن کاملRobust multigrid preconditioners for the high-contrast biharmonic plate equation
We study the high-contrast biharmonic plate equation with HCT and Morley discretizations. We construct a preconditioner that is robust with respect to contrast size and mesh size simultaneously based on the preconditioner proposed by Aksoylu et al. (2008, Comput. Vis. Sci. 11, pp. 319–331). By extending the devised singular perturbation analysis from linear finite element discretization to the ...
متن کاملMultilevel Preconditioners for Strongly Anisotropic Problems
In this dissertation, we develop new multilevel approaches to precondition algebraic problems stemming from the finite volume discretization of the diffusion equation with anisotropic, discontinuous coefficients. Two approaches are discussed. In the first approach, preconditioners are based on a partitioning of the mesh in the (x, y)-plane into non-overlapping subdomains and on a special coarse...
متن کاملMultigrid Preconditioners for Bi-cgstab for the Sparse-grid Solution of High-dimensional Anisotropic Diffusion Equation
Robust and efficient solution techniques are developed for high-dimensional parabolic partial differential equations (PDEs). Presented is a solver based on the Krylov subspace method Bi-CGSTAB preconditioned with d-multigrid. Developing the perfect multigrid method, as a stand-alone solver for a single problem discretized on a particular grid, often requires a lot of optimal tuning and expert i...
متن کاملA fast iterative solver for the variable coefficient diffusion equation on a disk
We present an efficient iterative method for solving the variable coefficient diffusion equation on a unit disk. The equation is written in polar coordinates and is discretized by the standard centered difference approximation under the grid arrangement by shifting half radial mesh away from the origin so that the coordinate singularity can be handled naturally without pole conditions. The resu...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Computat. and Visualiz. in Science
دوره 13 شماره
صفحات -
تاریخ انتشار 2010